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Abstract -Localization-delocalization transition in a discrete Anderson nonlinear Schrodinger 
equation with disorder is shown to be a critical phenomenon — similar to a percolation transition 
on a disordered lattice, with the nonlinearity parameter thought as the control parameter. In 
vicinity of the critical point the spreading of the wave field is subdiffusive in the limit t — > +co. 
The second moment grows with time as a powerlaw ex i°, with a exactly 1/3. This critical 
spreading finds its significance in some connection with the general problem of transport along 
separatrices of dynamical systems with many degrees of freedom and is mathematically related 
with a description in terms fractional derivative equations. Above the delocalization point, with 
the criticality effects stepping aside, we find that the transport is subdiffusive with a = 2/5 
consistently with the results from previous investigations. A threshold for unlimited spreading is 
calculated exactly by mapping the transport problem on a Cayley tree. 
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Introduction. — The problem of Anderson localiza- 
tion in disordered media is part of the general problem 
of transport of waves in disordered media. It came into 
focus after Anderson's suggestion that extensive disor- 
der traps electronic wave function in the tight binding 
regime [I]. In recent experiments with matter waves, ex- 
ponential localization of a Bose-Einstein condensate re- 
leased into a one-dimensional waveguide in the presence 
of a controlled disorder has been obtained [2J. The phe- 
nomenon is caused by destructive interference between the 
many scattering paths as the wave reflects off the struc- 
tural inhomogeneities of the medium. An interesting new 
topic came with the introduction of nonlinear models |3] 
and with the realization that a weak nonlinearity of the 
wave process can destroy localization (3J|4]. Indeed it was 
found in direct numerical simulations of one-dimensional 
discrete nonlinear Schrodinger lattice with disorder pjis] 
that above a certain critical strength of nonlinearity an 
unlimited spreading of the wave function occurs, with the 
second moment that grows with time as a powerlaw ex t" 
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and the exponent a in the range 0.3-0.4. The exact value 
of a has remained a matter of debate. 

In this Letter we analyze delocalization processes in dis- 
ordered media as a transport problem for a dynamical sys- 
tem with many coupled degrees of freedom. This problem 
is solved exactly on a Cayley tree. Our main findings are 
as follows. We show that the localization-delocalization 
transition is essentially a critical phenomenon — similar 
to a percolation transition in random lattices. Delocal- 
ization occurs spontaneously when the strength of nonlin- 
earity goes above a certain limit. We obtain the critical 
strength analytically and report it here for the first time. 
In vicinity of the critical point the spreading of the wave 
field is subdiffusive in the limit t — >■ -|-oo. The second 
moment grows with time as a powerlaw ex i", with a ex- 
actly 1/3. This critical regime is completely dominated by 
multi-scale long-time correlations, so that the correlation 
time is infinite. We describe this type of spreading in terms 
of a generalized diffusion equation with a fractional time 
the so-called Caputo fractional derivative. The critical 
regime finds its significance in connection with the general 
problem of transport along separatrices of dynamical sys- 
tems with many degrees of freedom |6j. Above criticality. 
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nonlinear properties take a stronger role over the dynam- 
ics. With the correlations effects stepping aside we find 
that the transport is subdiffusive with a — 2/5. We show 
that this regime corresponds to Markovian diffusion with 
a range-dependent diffusion coefficient consistently with 
the results from previous investigations (3l|4). This differ- 
entiation between the two transport regimes (1/3 versus 
2/5) is important at it helps to sort out some ambigui- 
ties in the reported transport exponents i4l[5] as well as 
to place the various transport models on a solid math- 
ematical background in connection with the asymptotic 
character of the transport. 

Model. — We work with a variant of discrete Ander- 
son nonlinear Schrodinger equation (DANSE) 



ih 



dtpn 
dt 



^LV'n+/3|V'„|Vr, 



with 



HL-ipn = Snipn + V'(V'n+l + V'n-l)- 



(1) 



(2) 



Here, H^ is the Hamiltonian of a linear problem in the 
tight binding approximation; /3 characterizes the strength 
of nonlinearity; on-site energies e„ are randomly dis- 
tributed with zero mean across a finite energy range; V 
is hopping matrix element; and the total probability is 
normalized to J^n IV'nP = 1- In what follows, ?i = 1 for 
simplicity. For /3 — > 0, DANSE (fTl) reduces to the origi- 
nal Anderson model Ij. All eigenstates are exponentially 
localized in that limit. We aim to understand the asymp- 
totic (t — > -f-cxi) spreading of initially localized wave packet 
under the action of nonlinear term. 

Analysis. — For each node n, we define a basis of 
linearly localized modes, the eigenfunctions of the lin- 
ear Hamiltonian Hl. By doing so we generate a map- 
ping of the wave number space into a functional space. 
For strong disorder, dimensionality of this space is infi- 
nite (countable). We consider this space as embedding 
space for the dynamicsr] If we introduce an inner prod- 
uct: tpn°''l'n = ^ni'^'nTi^n^ ^^ obtaiu a Hilbcrt space, in 
which an orthogonal basis can be defined. Here, the star 
denotes complex conjugate. Denoting the basis modes 
as (j)n,m, we have, with ujrn the eigenvalues of the linear 
problem, HL(j)n^m = ^m'Pn.m- Orthogonality implies that 
'}ln^*n,m4>n.k = Sm,k, whcrc Sm.k IS Kroncckcr's delta. In 
the basis of linearly localized modes the wave function ipn 
can be represented as 



tpn = /^0'm{t)4>n,-. 



(3) 



where we have introduced CTm(i), a set of complex am- 
plitudes which describe the evolution of the wave field 
in time. The total probability being equal to 1 implies 



^To visualize, think of each node as comprising a countable num- 
ber of compactified dimensions that are expanded to form a func- 
tional space. 



that J2m'^rnit)'^m{t) = 1. Wc now obtain a set of dy- 
namical equations for am{t). For this, substitute (l3| into 
DANSE (fll), then multiply the both sides by </»*_,„, and 
sum over n, remembering that the modes are orthogonal. 
The result reads 



-/3 Y. ^'=^"H. 



.CTm,, cr„,„cr,; 



ZCTfc — WfeCTfc — /J / ^ i'/c,mi,m2,m3"mii^„i2"'"3' 

nil ^^n2^m,3 



(4) 



where the coefficients Vk.rni,m2.m3 are given by 

<'k,mi,m2,m3 — / ^ yn.fcyw.Tti yrt.m^ y»,™3 ' \^) 

n 

and we have used dot to denote time derivative. Equa- 
tions (Hi) correspond to a system of coupled nonlinear os- 
cillators with the Hamiltonian 



H = Hn + Hi, 



Ho-T. 



^k<Jk^k, 



(6) 



int — r, /_^ Vk,mi.m2,m3<^k'^rni<^m2'^rna- (') 

k^mi ,m2,rM3 

Here, Hq is the Hamiltonian of non-interacting harmonic 
oscillators and Hint is the interaction Hamiltonian|j Each 
nonlinear oscillator with the Hamiltonian 

hk = ^k<y*k^k + -^yk,k:kM<^*k^k(^*k^k (8) 

and the equation of motion 

idk - uJkCTk - PVk,k,k,kCrkCrlcrk = (9) 

represents one nonlinear eigenstate in the system — iden- 
tified by its wave number fc, unperturbed frequency Wfc, 
and nonlinear frequency shift Acj^ = f3Vk^k.k,kO'kO'k- Non- 
diagonal elements Vk^mi,m2.m3 characterize couplings be- 
tween each four eigenstates with wave numbers k, mi, TO2, 
and m^. It is understood that the excitation of each eigen- 
state is not other than the spreading of the wave field in 
wave number space. Resonances occur between the eigen- 
frequencies Uk and the frequencies posed by the nonlinear 
interaction terms. We havejj 



Wfe 



(10) 



When the resonances happen to overlap, a phase trajec- 
tory may occasionally switch from one resonance to an- 
other. As Chirikov realized 7|, any overlap of resonances 
will introduce a random element to the dynamics along 
with some transport in phase space. Applying this argu- 
ment to DANSE (fll) , one sees that destruction of Anderson 
localization is limited to a set of resonances in a Hamil- 
tonian system of coupled nonlinear oscillators, eqs. (pi) 
and ([7|), permitting a connected escape path to infinity. 

^We include self-interactions into -ffint- 

^Conditions for nonlinear resonance are obtained by accounting 
for the nonlinear frequency shift. 
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At this point, the focus is on topology of the random 
motions in phase space. We address an ideahzed situa- 
tion first where the overlapping resonances densely fill the 
phase space. This is fully developed chaos, a regime that 
has been widely studied and discussed in the literature 
(e.g., Refs. |8Jl9]). Concerns raised over this regime when 
applied to eqs. ([6| and ([t]) come from the fact that it re- 
quires a diverging free energy reservoir in systems with a 
large number of interacting degrees of freedom. Yet, devel- 
oped chaos offers a simple toy-model for the transport as 
it corresponds with a well-understood, diffusive behavior. 

A more general, as well as more intricate, situation oc- 
curs when the random motions coexist along with regular 
(KAM regime) dynamics. If one takes this idea to its ex- 
treme, one ends up with the general problem of transport 
along separatrices of dynamical systems. This problem 
constitutes a fascinating nonlinear problem that has as 
much appeal to the mathematician as to the physicist. 
An original important promotion of this problem to large 
systems is due to Chirikov and Vecheslavov [6]. 

This type of problem occurs for slow frequencies. One 
finds 



10 11 



that resonance-overlap conditions are satis- 
fied along the "percolating" orbits or separatrices of the 
random potential where the orbital periods diverge. The 
available phase space for the random dynamics can be 
very "narrow" in that case. In large systems, the set of 
separatrices can moreover be geometrically very complex 
and strongly shaped. Often it can be envisaged as a frac- 
tal network at percolation as for instance in random fields 
with sign-symmetry [12| . 

There is a fundamental difference between the above two 
transport regimes (chaotic versus near-separatrix) . The 
former regime is associated with an exponential loss of 
correlation permitting a Fokker-Planck description in the 
limit t — ?> -fcxD. The latter regime when considered for 
large systems is associated with an algebraic loss of corre- 
lation, implying that the correlation time is infinite. There 
is no a conventional Fokker-Planck equation here, unless 
generalized to fractional derivatives [l3{[l4] , nor the famil- 
iar Markovian property (i.e., that the dynamics are mem- 
oryless). On the contrary, there is an interesting inter- 
play fll] between randomness, fractality, and correlation, 
which is manifest in the fact that all Lyapunov exponents 
vanish in the thermodynamic limit, despite that the dy- 
namics are intrinsically random. 

This situation of random non-chaotic dynamics with 
zero Lyapunov exponents, being in fact very general 
15 16 , has come to be known as "pseudochaos." One 
might think of pseudochaos as occurring "at the edge" 
of stochasticity and chaos, thus separating fully developed 
chaos from domains with regular motions. In what follows, 
we discuss the implications of chaotic and pseudochaotic 
transport for the spreading of wave field in eqs. (4). 

Chaotic case. — As the time correlations vanish expo- 
nentially fast, equation (4) can conveniently be considered 
as a Langevin equation with the nonlinear interaction term 
thought as a white noise term in the limit t — >■ -|-c». There 



is a well-defined diffusion coefficient here, which behaves 
as modulus square of the complex interaction amplitude. 
The cubic interaction in eq. B implies that D c>c |cr„|^. 
If the field is spread over An sites, then the conserva- 
tion of the probability dictates |cr„P ~ 1/An, leading to 
D ex 1/{Ati)^. One sees that the transport problem in the 
chaotic case is essentially a diffusion problem with a range- 
dependent diffusion coefficient. This range-dependence is 
an inverse powerlaw, as one would indeed expect for a 
homogeneous random system. 

Let / — f{t, An) be the probability density to find an 
initially localized wave field at time t at distance An from 
the origin. The diffusive character of the spreading justi- 
fies 



d fit, An) 
dt 



d 



dAn 



W 



1 df{t,An) 
(An)3 dAn 



(11) 



where IF is a constant coefficient which collects all rele- 
vant parameters of the diffusion process. The fundamental 



solution or Green's function of eq. (11 1 reads 



f{t,An) 



r(l/5) (25VFt)i/5 



exp 



(An)5 
2bWt 



(12) 



where the normalization J„ f{t, An)dA7i — 1 has been 
applied. It is noticed that the distribution in eq. ( 12 ) is 
essentially non-Gaussian. From eq. ( |12[ ) one immediately 
obtains that 

{{An)^(t)) = (25m)2/5r(3/5)/r(l/5), (13) 

where the angle brackets denote ensemble average, and F is 
the Euler gamma function. The net result is ((An)^(i)) ex 
i^/^, consistently with the finding of Refs. (slkj. 

Pseudochaotic case. — This regimes takes eqs. Q to 
an opposite extreme limit where each oscillator can only 
communicate with the rest of the wave field via a nearest- 
neighbor rule. This is a marginal regime yet permitting an 
escape path to infinity. Glearly, the number of coupling 
links is minimized in that case. When summing on the 
right-hand-side, the only combinations to be kept are, for 
the reasons of symmetry, crfccr^crfe and ak^ia^ak+i. We 
have 

i&k - LOkCTk = PVkO-kcrlak + 2(iV^ak-i'yl(Tk+i, (14) 



where we have also denoted for si mpli city Vk = Vk 



k^k,k 

and V^ = Vk^k-i,k.k+i- Equations (14) define an infinite 



(fc = 1, 2, . . .) chain of coupled nonlinear oscillators where 
all couplings are local (nearest- neighbor- like) . The inter- 
action Hamiltonian in cq. ([7]) is simplified to 



Hu 



/3 



E^^- 



kO-k(Tk(Tk ■ 



l^22^k''^*k<^k-lcr*k<Jk+l- 



(15) 
We are now in position to introduce a simple lattice model 



for the transport. The key step is to observe that eqs. ( 14 ) 
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Fig. 1: Mapping eqs. (14 1 on a Cayley tree. Each node rep- 



resents a nonlinear eigenstate, or nonlinear oscillator with the 
equation of motion idfc — ijJ-k(yh — l3Vk,k,k,kO'kO'k<^k = 0. Blue 
nodes represent oscillators in a chaotic ("dephased") state. 
Black nodes represent oscillators in regular state. One ingo- 
ing and two outgoing branches on node k {k — 1,2, . . .) repre- 
sent respectively the complex amplitudes a^, ak-i, and ak+i- 
Structures that are not explicitly shown are beyond the dashed 
lines. 



can be mapped on a Cayley tree where each node is con- 
nected to 2 = 3 neighbors (here, z is the coordination 
number). The mapping is defined as follows. A node with 
coordinate k represents a nonlinear eigenstate, or nonlin- 
ear oscillator with the equation of motion ^ . There are 
2 = 3 branches at each node: one that we consider ingoing 
represents the complex amplitude ct^, and the other two, 
the outgoing branches, represent the complex amplitudes 
CTfe-i and CTfc+i respectively. These settings are schemati- 
cally illustrated in Fig. 1. 

A Cayley tree being by its definition 
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cal graph offers a suitable geometric model for infinite- 
dimensional spaces. We think of this graph as embedded 
into phase space of the Hamiltonian system of coupled 
nonlinear oscillators, eqs. ^ and (l7|. In the thermody- 
namic limit fcmax — >■ oo, in place of a Cayley tree, one uses 
the notion of a Bethe lattice]^ Setting /cmax — >■ oo, we sup- 
pose that each node of the Bethe lattice hosts a nonlinear 
oscillator, eq. ^. The bonds of the lattice, in their turn, 
can conduct oscillatory processes to their neighbors as a 
result of the interactions present. 

Next, we assume that each oscillator can be in a chaotic 
( "dephased" ) state with the probability p (and hence, in a 
regular state with the probability 1 — p). The p value 
being smaller than 1 implies that the domains of ran- 
dom motions occupy only a fraction of the lattice nodes. 
Whether an oscillator is dephased is decided by Chirikov's 
resonance-overlap condition — which may or may not be 



*A Bethe lattice is an infinite version of the Cayley tree. To this 
end, a purist might prefer to say "bond" in place of "branch," but 
that's all about the terminology. 



matched on node k. We believe that in systems with many 
coupled degrees of freedom each such "decision" is essen- 
tially a matter of the probability. The choice is random. 
Focusing on the p value, we consider system-average non- 
linear frequency shift 



AWNL = /^dV"?! 



/An 



(16) 



as an effective "temperature" of nonlinear interaction. It 
is this "temperature" that rules over the excitation of the 
various resonant "levels" in the system. With this inter- 
pretation in mind, we write p as the Boltzmann factor 



p — exp(— (Jw/Aojnl) 



(17) 



where 5uj is the characteristic energy gap between the res- 
onances. Expanding ipn over the basis of linearly localized 
modes, we have 



(1^. 



A 



v„ Z.^ l-^ 



(18) 



n nil ,ni2 



The summation here is performed with the use of orthog- 



onality of the basis modes. Combining with eq. (16), 



Awnl 



All 



Y^ 



(7,y,0-r, 



/ ,"m^™ 



(19) 



The sum over m is easily seen to be equal to 1 due to 
the conservation of the probability. Thus, Awnl — /3/An. 
If the field is spread over An states, then the distance 
between the resonant frequencies behaves as Stu '^ l/An. 
We normalize units in eq. (fTl) to have Stu — l/An exactly. 
One sees that p — exp(— 1//3). This result suggests that 
behavior be non-perturbative in the pseudochaotic regime. 
For the vanishing /3 — > 0, the Boltzmann factor p — > 0, 
implying that all oscillators are in regular state. In the 
opposite regime of /? — ?> cxd, p — > 1. That means that 
all oscillators are dephased and that the random motions 
span the entire lattice. 

There is a critical concentration, pc, of dephased oscil- 
lators permitting an escape path to infinity for the first 
time. This critical concentration is not other than the 
percolation threshold on a Cayley tree. In the basic the- 
ory of percolation it is found that p^ — l/(z — 1) (e.g., 
Ref. [it]). This is an exact result. For z = 3, Pc = 1/2. 
We associate the critical value pc = 1/2 with the onset of 
transport in the DANSE model, eq. (II]). When translated 
into the f3 values the threshold condition reads 



/3c = l/ln(z-l). 



(20) 



Setting z = 3, we have /3c = l/ln2 w 1.4427. This value 
defines the critical strength of nonlinearity that destroys 
the Anderson localization. For the /3 values smaller than 
this, the localization persists, despite that the problem is 
nonlinear. When /? > l/ln2, the localization is lost, and 
the wave field spreads to infinity. 
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Our conclusion so far is that the loss of localization is 
a threshold phenomenon, which requires the strength of 
nonlinearity be above a certain level. In this respect, the 
nonlincarity parameter /3 acquires the role of the control 
parameter. The onset of unlimited spreading is at /3c = 
l/ln2, this being an exact result of the model. We now 
turn to predict second moments for the onset spreading. 

This task is essentially simplified if one visualizes the 
transport as a random walk over a system of dephased 
oscillators. For p ^ pc, this system is self-similar, i.e., 
fractal. That means that dephased oscillators form ar- 
bitrarily large clusters, each presenting the same fractal 
geometry of the infinite percolation cluster 18 . It is the 



infinite cluster that conducts unlimited spreading of the 
wave function on a Bethe lattice. We should stress that 
the fractal geometry of the clusters is a consequence of the 
probabilistic character of dephasing. 

In random walks on percolation systems one writes the 
mean-square displacement from the origin as 
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((An)2(i))=A2i2/(2+e)^ t^+oo. 



(21) 



Here, ^ is a dimensional constant parameter and 6 is the 
exponent of anomalous diffusion, or the connectivity ex- 
ponent. This exponent accounts for the deviation from 
the usual Fickian diffusion in fractal geometry. As a rule, 
9 > 0, implying that the diffusion is slowed down in com- 
plex "labyrinths" of the fractal. Another common way of 



writing eq. (21 ) is given by 18 



((An)2(t))=^2^'^,/d,^ i^+oo. 



(22) 



where df is the Hausdorff dimension which measures the 
number of nodes that belong to a given cluster, and ds = 
2df /{2 + 0) is the fracton, or spectral, dimension which 
describes the density of states in fractal geometry. The 
key difference between the two dimensions is that df is a 
purely structural characteristic of the fractal, whereas ds 
reflects the dynamical properties such as wave excitation, 
diffusion, etc. Note, also, that the spectral dimension is 
not larger than df. 

The two scaling laws above apply to any percolation 
system. For percolation on a Cayley tree, the following 
exact results hold 



One sees that 



18 
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= 4, d/ = 4, and 4 = 4/3. 



((An)2(i)) cxt^/^ t 



-l-oo. 



(23) 



This is the desired scaling. By its derivation, subdiffusion 
in eq. (231 is asymptotic (t — >■ -foo) in the thermodynamic 
limit. 

We proceed with a remark that the Hausdorff dimension 
being equal to 4 matches with the implication of eqs. m 
and ([5]) where the coefficients Vk.mi,m2.m3 are supposed 
to run over 4-dimensional subsets of the ambient Hilbert 
space. Indeed it is the overlap integral of four Ander- 
son eigenmodes, eq. ([S]), that decides on dimensionality of 
subsets of phase space where the transport processes con- 
centrate. When the nearest-neighbor rule is applied, this 



overlap structure is singled out for the dynamics. Under 
the condition that the structure is critical, i.e., "at the 
edge" of permitting a path to infinity, the support for the 
transport is reduced to a percolation cluster on a Bethe 
lattice — characterized, along with the above value of the 
Hausdorff dimension, by the very specific connectivity ex- 
ponent, = 4. The end result is a = 2/(2 + 9) = 1/3. 



At contrast with eq. (11 1, random walks on percolation 



systems are described by a non-Markovian diffusion equa- 
tion |11[|21|: 



d fit, An) 
dt 



df 



d d^ fit', An) 



re Jo it - t')<^/i^+'^) df 9(An)2 



(24) 



where Tg = r(2/(2 -I- 9)). In writing eq. ([24]) we have 
adopted results of Ref. [Til to random walks on a single 
cluster. The convolution on the right hand side is express- 
ible in the compact form of a fractional time the so-called 
Caputo fractional derivative 22 of order < 9/{2+9) < 1. 



We have, in the commonly used notations. 



d fit. An) 
dt 



9/(2+9) 9^ fit. An) 
° * a(An)2 



(25) 



Allowing for 9 = 4, one sees that the critical spreading 
corresponds to a fractional operator qD^. = qDI~" for 
a = 1/3. 

The latter result should be addressed. Indeed eq. (251 



or "strange,' 



shows that the critical spreading is a matter of fractional, 
kinetics 23 , consistently with the implica- 

In many ways 
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tion of pseudochaotic behavior 
equations built on fractional derivatives offer an elegant 
and powerful tool to describe anomalous transport in com- 
plex systems. There is an insightful connection with a 
generalized master equation formalism along with a math- 
ematically convenient way for calculating transport mo- 
ments as well as solving initial and boundary value prob- 
lems 14 24 . The fundamental solution of the fractional 
eq. (25 1 is evidenced in Table 1 of Ref. [24 . It shares non- 



Gaussianity with Green's function in eq. ( 12 ), being in the 
rest analytically very different. 

Recently, an application of the fractional diffusion equa- 
tion to subdiffusion in the nonlinear Schrodinger equation 
with disorder has been proposed within a continuous time 
random walk (CTRW) formalism 25 . It was shown that 



there is a heavy-tailed distribution of waiting times be- 
tween consecutive steps of the random motion, thus giv- 
ing rise to a slower-then-linear growth of second moments. 
In the above eqs. (21) and (23) we have not as a matter 



of fact assumed any heavy-tailed distribution of this sort. 
At contrast, in the present model, the random walker is 
supposed to take one unit step along the cluster as soon 
as one unit time is elapsed. Despite this difference, the 
two approaches are essentially equivalent and do lead to 
the same type of fractional diffusion equation, eq. (25). 



This is because fractal labyrinths of the subset that hosts 
the random dynamics act as to introduce multi-scale trap- 
pings to the motion characterized by an inverse powerlaw- 
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like distribution of waiting times. The latter is found to 
be, in properly normalized units, 

0(r)ocl/(l + r)(4+^)/(2+«). (26) 

As mean effective waiting time diverges, the basic assump- 
tions of CTRW's are readily installed. 

Summary. — We have shown that the Anderson lo- 
calization in disordered media can be lost in the presence 
of a weak nonlinearity and that the phenomenon is critical 
(thresholded). That means that there is a critical strength 
of nonlinearity above which the wave field turns to an un- 
limited spreading. Below that limit, the field is localized 
similarly to the linear case. We have discussed the problem 
as a percolation problem on a separatrix system of discrete 
nonlinear Schrodinger equation with disorder. This prob- 
lem is solved exactly on a Bethe lattice. A threshold for 
delocalization is found to be /3c = l/ln2 « 1.4427. 

For the /3 values smaller than this, the localization per- 
sists, despite that the problem is nonlinear. Support for 
this type of behavior can be found in the results of Ref. 



26 . In vicinity of the delocalization point the spread- 
ing of the wave field is subdiffusive, with second moments 
that grow with time as a powerlaw ex t" for t — > +oo. 
We find that a is exactly 1/3 in the thermodynamic limit. 
This regime bears signatures enabling to associate it with 



the onset of "weak" transport 27 of Alfven eigenmodes 



in vicinity of marginal stability of magnetic confinement 
systems. 

The key feature of the critical spreading is that it is com- 
pletely dominated by multi-scale time correlations that 
are long-ranged. We associate this regime with Hamil- 
tonian pseudochaos, random non-chaotic dynamics with 
zero Lyapunov exponents 13 15 16 . The implications of 



pseudochaos find their significance in the general picture 
of transport along separatrices of Hamiltonian systems, 
the cornerstone of nonlinear dynamics, and are mathe- 
matically related with a description in terms of fractional 
derivative equations. We find it interesting to remark 
that the above fractional diffusion equation, eq. (25), is 
"born" within the mathematical structure of nonlinear 
Schrodinger equation with usual time differentiation. It is 
in fact the interplay of nonlinearity and randomness that 
leads to a fractional derivative equation when it comes to 
non-Markovian transport of the wave field. 

Above criticality, as /3 increases, nonlinearity takes a 
stronger role over the dynamics. To this end, the multi- 
scale and pseudochaotic properties turn to lose their im- 
portance. For the large /3 ^ 1, one may go for a toy- 
model that singles out the nonlinear properties of the 
transport. One such model is a Markovian diffusion model 
with range-dependent diffusion coefficient. In the latter 
case we find a = 2/5 consistently with the results from 
previous investigations (slli]. 
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